Collective excitations of trapped one-dimensional dipolar quantum gases 
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We calculate the excitation modes of a ID dipolar quantum gas confined in a harmonic trap with 
frequency u>o and predict how the frequency of the breathing n = 2 mode characterizes the interaction 
strength evolving from the Tonks-Girardeau value lo 2 = 2lo to the quasi-ordered, super-strongly 
interacting value ui2 = V^loo- Our predictions are obtained within a hydrodynamic Luttinger-Liquid 
theory after applying the Local Density Approximation to the equation of state for the homogeneous 
dipolar gas, which are in turn determined from Reptation Quantum Monte Carlo simulations. They 
are shown to be in quite accurate agreement with the results of a sum-rule approach. These effects 
can be observed in current experiments, revealing the Luttinger-liquid nature of ID dipolar Bose 
gases. 
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Introduction.- The bottom line in the design of con- 
ceptually novel technological applications is the possi- 
bility of reaching extreme quantum degeneracy under 
controlled conditions. This is quite a remarkable prop- 
erty of ultracold atomic gases, which is being evidenced 
by the interdisciplinary contribution of quantum optics, 
quantum information, condensed matter, atomic, and 
fundamental physics. Extreme quantum limit can here 
be obtained after lowering the temperature down to the 
nanokelvin scale, and by tuning the atomic interactions 
and the dimensionality. Dimensionality can indeed be re- 
duced to two and one dimensions (ID) after using a vari- 
ety of magnetic and optical techniques including optical 
lattices [ITS 0, 0] • The interactions can also be manipu- 
lated almost at will in their short-range part by means of 
the Fano-Feshbach resonance mechanism [a, la, w, la la] • 
In the case of atomic or molecular species with large mag- 
netic or electric moments, proposals have been put for- 
ward which predict the possibility of tu ning the long- 
range dipolar tail of the interaction 0, \lM- Th e ob- 
servation of dipolar interactions [13] in atomic 52 Cr va- 
pors with relatively large magnetic moment \x ~ 6jUs, 
is especially promising for applications, since it opens 
the way to the realization of e.g. different quantum 
phases I13l |l4j and the observation of spin-charge sep- 



aration [15[. Molecular dipolar crystals have been more 
recently proposed as a realization of high fidelity quan- 
tum memory for quantum computation [161 ]. In a set of 
landmark experiments 13, HI and supported by simu- 
lational work [l9l_|20|, vapors of 52 Cr atoms have been 
Bose-condensed [21]. Yet, the effect of the dipolar in- 
teraction has been largely enhanced after reducing the 
strength of the short-range part [22| • 

Combination of ID geometries with tunable interac- 
tions may provide an easy access to enhanced quantum 



correlations. Experimental realizations of the strongly 
correlated Tonks-Girardeau (TG) gas [23| have already 
been observed in atomic vapors with contact interac- 
tions [2J, [25[ . Going beyond the TG regime under these 
conditions is not an easy task (2(| . However, as we have 
more recently demonstrated by a combined Reptation 
Quantum Monte Carlo (RQMC) and bosonization ap- 
proach [27(, a homogeneous ID Bose gas with dipolar 
interactions is a very strongly correlated Luttinger liquid 
with the parameter K < 1 [28( at all densities nro ^ 0.1 
with ro being the range of the dipolar potential (see be- 
low). The Luttinger liquid crosses over from a Tonks- 
Girardeau gas setting in for nro 0-1 to a high-density 
quasi-ordered state [2 9j which can be viewed as the ana- 
logue of a Charge Density Wave. We have also discussed 
how the use of polar molecules may provide access to 
this quasi-ordered state, which from now on we refer to 
as Dipolar Density Wave (DDW). 

The knowledge of the interaction regime is a basic tool 
for all conceivable applications. Since ongoing experi- 
ments will be performed in confined geometry, one of the 
best suited and controlled methods is to exploit the col- 
lective excitations related to the discrete modes in the 
harmonic trap [301 ] . which has indeed been largely used 
since the very first experiments 31 , 3^] • The study of the 
collective modes can also be useful in view of the proposal 
to implement the quantum memory, to reveal and inves- 
tigate decoherence mechanisms possibly arising from the 
coupling of internal and external degrees of freedom on 
the DDW state [3]. 

In this Letter we predict for the first time the crossover 
behavior of the collective modes of the trapped ID dipo- 
lar Bose gas by two different theoretical methods, namely 
a sum-rule approach and a hydrodynamic Luttinger- 
liquid model. The two methods share the application 



2 



of the Local Density Approximation to the equation of 
state as determined from our RQMC simulational data, 
and are shown to give results in quite remarkable agree- 
ment. 

The equation of state.- We first determine the 
ground-state energy per particle (within the statisti- 
cal error) of the homegeneous dipolar Bose gas by re- 
sorting to Reptation Quantum Monte Carlo simulations 
(RQMC) (33)]. As described in more detail in Ref. [13], 
we consider TV atoms or molecules of mass M and 
permanent dipole moments arranged along a line in 
the limit of negligible contact interaction [22| and po- 
larized in the orthogonal direction. The Hamiltonian 

H = (-(nr ) 2 E(3 2 /^ 2 ) + Ko) 3 E l<J l^-^l 3 ) ^ 
defined in effective Rydberg units Ry* = H 2 / (2Mr 2 ) with 
r = MCdd/ (2ttH 2 ) the effective Bohr radius, Cdd — MoM 2 
or Cdd = d 2 /eq the interaction stengths for magnetic fid 
or electric d dipole moments respectively. The governing 
dimensionless parameter is nro, with n the number of 
dipoles per unit length. 

The RQMC data for the dependence of the energy per 
particle e on n are fitted by the expression e(nro)/ Ry* = 
[C(3)(nr ) 4 + a(nr ) e + b(nr ) f + c(nr ) 2+9 } [1 + Tiro]" 1 + 
(7r 2 /3)(nr ) 2 [l+d(?ir ) 9 ]- 1 , where a = 3.1(1), b = 3.2(2), 
c = 4.3(4), d = 1.7(1), e = 3.503(4), / = 3.05(5), 
and g — 0.34(4) with an overall x 2 ed — 5- For 
nr Q -C 1, e(nr )/Ry* ~ (it 2 /3)(w ) 2 and for nr > 1, 
e{nr )/Ry* ~ C(3)(rw ) 3 so that both the TG and DDW 
limits are satisfied. The functional form for e(n) here 
provided can be used in further calculations. Here, we 
use it to obtain the chemical potential /^rqmc(^) = 
(l + n(d/dn))e(n). 

The Local Density Approximation.- In experiments, 
the effectively ID dipolar gas is confined by a harmonic 
potential V(x) — Moj 2 x 2 /2 in the axial direction. The 
dynamical behavior of the confined gas can be found from 
hydrodynamic equations at temperature T = using a 
Local Density Approximation (LDA) to the equation of 
state, in which the energy of the inhomogeneous system 
is a local functional of the local density n(x) expressed 
as the integral over x of the energy density of the homo- 
geneous gas with density n(x). The validity of such an 
approach is limited to dynamical effects in which the typ- 
ical length over which n(x) varies is much larger than the 
average interparticle distance in the axial direction. Fur- 
thermore, the dynamical behavior in the radial (trans- 
verse) plane must remain frozen. Under these conditions, 
the ground-state density profile n(x) of the trapped dipo- 
lar gas is obtained by plugging /xrqmc(^) into the equa- 
tion 



H{n(x)) + V(x) = no , 



(1) 



and solving it for fig and n(x) with the condition that the 
total number of particles N = J^ R n(x)dx is conserved, 
where the Thomas-Fermi radius R = (2ho/(Muj 2 )) 1 / 2 is 



such that n(±i?) = 0. The condition can be casted in the 
form N(r /a ho ) 2 = ( Mo /Ry*) 1 ' 2 /i r^- 1 ^/ Ry*){\ - 
x 2 )]dx, where a^ = (H/Mlliq) 1 ^ 2 is the harmonic oscilla- 
tor length. This identifies N(ro/aho) 2 as the interaction 
parameter driving the trapped dipolar gas from the TG 
across the DDW regime. 

The evolution of the calculated density profiles through 
the crossover, shows the expected increase of the central 
density n(0) with the chemical potential /xo, as well as a 
steepening of the profiles at the trap edges. The calcu- 
lated profiles agree with the analytical results expected 
in the TG and in the DDW limits. 

The breathing mode from a sum rule (SR) approach.- 
We first determine the effect of the interactions on 
the frequency of the lowest compressional (breathing) 
mode of the trapped gas by a sum-rules approach. This 
mode is coupled to the ground state by the operator 
X2 = ^2iLi x f- This makes it the easiest to probe 
in current experiments, as it is excited by modulat- 
ing the trap frequency and observed by following the 
time-evolution of the width of the cloud by conventional 
density- imaging techniques. The frequency of the 
breathing mode satisfies an inequ ality that can be de- 
rived from the sum rule [3, 13a , l3(j ] us B < mi/m_i, 
where 2mi = ([X, [H, X]}) and 2m_i is the static re- 
sponse function. After some simplification, this yields 
the upper bound: 



< ni 



d{x 2 )/duj 2 



(2) 



where (x 2 ) — N^ 1 J x 2 n(x)dx. A closed form expression 
of VLb in Eq. j2]) can be obtained by means of a scaling 
argument whenever fi(n) is of the form /i(n) = An 7 , re- 
sulting in £Ib = ujo(2 + 7) 1 / 2 . In particular, in the TG 
case with 7 = 2 this yields Clg = 2ujq, and in the DDW 
case with 7 = 3 it yields f2s < V^o- 

For intermediate interaction strengths, we have to re- 
sort to a numerical estimation of Eq. @ using the LDA 
density profile. The result is represented by the solid line 
in Fig. [l] showing the smooth evolution of the breathing 
mode frequency from the TG to the DDW regimes. 
Figure [1] represents an useful tool to identify the inter- 
action regime of the trapped dipolar gas through one 
of the best handled experimental probes available with 
cold atomic (molecular) quantum gases 3l|, [32j. How- 
ever, Eq. ([2]) in principle yields an upper bound for the 
frequency of the lowest compressional mode. Moreover, 
estimating the frequencies of the higher modes by the 
sum-rule approach is not as simple. We thus switch to 
an alternative, hydrodynamic, approach to compute their 
frequencies. 

The excitation modes from hydrodynamic Luttinger 
equations.- In [271 ]. we have shown that the low-energy 
behavior of a homogeneous dipolar Bose gas is well de- 
scribed by the Luttinger hamiltonian, and we have de- 
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FIG. 1: Squared frequency u>% of the breathing mode scaled 
to the trap frequency cjq vs. the interaction parameter 
N(ro/a,ho) 2 , as calculated from two models. Solid line: sum- 
rule approach Eq. @ (SR). Symbols: Luttinger-liquid hydro- 
dynamics Eq. ((6]) (LL). 



termined the density-dependence of the velocity u and 
Luttinger exponent K by combining bosonization and 
RQMC techniques. The u and K obtained from the 
RQMC structure factor were found to agree with those 
extracted from the RQMC energy via the Luttinger rela- 
tions uK = M _1 7T/ 



Galilean invariance 



= ■K~ l d n ^{n) embodying 
We now assume that in a 
slowly varying external trapping potential, the dipolar 
Bose gas can be described by a Luttinger liquid hamilto- 
nian, 



H 



LL 



dx 
2^ 



u(x)K(x){ttU) 2 



u{x) 
K{x) 



(3) 



where u{x) and K(x) now depend on position via the 
LDA n(x). They are related through u(x)K(x) = 
M irn(x) and u(x)/K(x) = w d n fi(n)\ n=n r x \, extend- 
ing the Luttinger-liquid relations to the weakly inhomo- 
geneous system. In Eq. ((3]), <j> and II satisfy canonical 
commutation relations [0(x),IT(y)] = iS(x — y), leading 
to the equations of motion: 



d t <j>(x,t) 
Trd t U(x,t) = d x 



ttu(x)K(x)TI(x, t), 
u(x) 



K(x) 



d x cj>(x,t) 



(4) 
(5) 



These equations of motion must be complemented by 
boundary conditions expressing that no current flows 
across the edges, i. e. j(±R, t) = 0. Since Sn = —d x (f>/Tr, 
the continuity equation dt (Sn) + d x j = leads naturally 
to j = dt4>/^, and allows us to rewrite the boundary 
conditions as <ft{—R) = 0o and <fi(R) = </>i (3^ |. 

Eqs. (HI)-© possess a stationary solution, d x 4> oc 
K(x)/u(x). This solution actually describes the addition 



of one particle to the system cj)(x) — > (j)(x) — -Kg(x)/g(R) 
where g(x) = J_ R (K(y)/u(y))dy. A consequence is that 
</>o and (j)i are not independent, but related through 
4>i — (fio = —ttN with N being the number of parti- 
cles added to the system. This form generalizes the 
bosonization formula derived in the case of a homoge- 
neous system [37]]. Alternatively, the static solution can 
be derived by considering the effect of a perturbation 
to the trapping potential [ifj| . Combining Eqs. (j!])-© 
and using linearity to search for solutions of the form 
4>(x, t)=(j> - TrNg{x)/g{R) + £„ >x Vn {x)e^ 1 we find: 



Muj lVn = n(x)d x (d n /j,(n(x))d x tp n ) , 



(6) 



with the boundary conditions tp n (±R) — for the dis- 
crete Fourier components ip n . 

While Eq. © is cast in a form identical to the hydro- 
dynamic equation for density excitations (3f3 | , we remark 
that here it has been derived from Hll Eq. |3j). Thus, a 
comparison of the measured excitation frequencies with 
those predicted from ^ with our RQMC data, provides 
a test of Luttinger-liquid behavior in trapped ID dipolar 
Bose gases within the validity of LDA. 

In the case of a harmonic trapping, one of the eigen- 
frequencies in Eq. © is obtained straightfowardly. In- 
deed, substituting fi(x) oc n{x) in Eq. ©, differentiat- 
ing Eq. © with respect to x and using the Luttinger- 
liquid relations, we find that ^>\{x) is an eigenfunction 
of © associated with the eigenvalue ujq. This particu- 
lar solution is simply the Kohn mode or sloshing mode 
describing a center-of-mass oscillation. Indeed, it can be 
recovered by expanding to first order in A the expression 
4>{x,t) = <fi(x — A cos(u>ot), 0) which describes a center of 
mass motion with rigid density. The eigenfrequencies are 
exactly known also in the two asymptotic limits. Indeed, 
insertion of oc n 1 in Eqs. © and ©, yields solutions 
of the form <p n (x) = A n (l - x 2 /R 2 fhc£ h+1/2) (x/R), 
the associated eigenvalues being cu 2 = luq7i[2+ (n— 1)7]/2 
, where Cn ^ a re G egenbauer polynomials and A n 
normalization factors [4l[, Thus, at the two opposite TG 
(7 = 2) and DDW (7 = 3) limits one finds respectively 
uj n = n 2 ujQ and ui n — n(3n — 1)o->q/2. For intermedi- 
ate densities, we have solved Eq. © using the Sledge 
algorithm available online [42| , after inserting as ingre- 
dients the computed LDA density profiles from Eq. © 
and of the analytical expression for d n fi obtained from 
the RQMC fit, evaluated at the local density for differ- 
ent values of the interaction parameter N(ro/aho) 2 - In 
the numerical solution we have taken special care of the 
finite mesh-size effects for best accuracy. 

The values of the breathing frequency wg/wo obtained 
from Eq. © are represented in Fig. Q] by the filled sym- 
bols, and agree up to the second digit with the sum- 
rule result. This was expected in the TG and DDW 
regimes where, as already noticed in Ref. 36], an equal- 
ity sign holds in Eq. ©. The eigenfrequencies of the 
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FIG. 2: (u) n /Lj ) 2 vs. N(r /a aho ) 2 from Eq. ©. From bottom 
to top: the modes with n = 2, 3, 4, and 5. Tick solid lines: 
limiting values in the TG and DDW regimes at each n. 



higher modes with n = 3, 4, and 5 are displayed in Fig. [5] 
as functions of N(ro/aho) 2 , showing the same smooth 
crossover behavior between the two opposite TG and 
DDW regimes. The exact frequencies in these asymp- 
totic regimes are recovered by the numerical calculations. 

Conclusions.- In conclusion, we have predicted the 
evolution of the collective modes of a ID dipolar Bose gas 
through the crossover from the Tonks- Girardeau to the 
Dipolar-Density- Wave regime. These modes, and espe- 
cially the breathing mode can be excited and measured 
by quite standard and reliable techniques. Our results 
are relevant to experiments, where they can be used to 
determine the interaction regime, and extended to inves- 
tigate the occurrence of decoherence effects in quantum 
applications . From the theoretical point of view, our 
results confirm for the trapped gas the smooth crossover 
from the TG to the DDW regimes 



27| . as expected 



for weak inhomogeneity. The RQMC functional form of 
the energy per particle that we explicitly provide in this 
work, can be useful to investigate the issue, crucial for 
applications, of the crystal phase stabilization by means 
of e.g. a commensurate, though shallow, optical lattice. 
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